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Abstract 

The convergence properties of numerical Regge calculus as an approximation to con¬ 
tinuum vacuum General Relativity is studied, both analytically and numerically. The 
Regge equations are evaluated on continuum spacetimes by assigning squared geodesic 
distances in the continuum manifold to the squared edge lengths in the simplicial man¬ 
ifold. It is found analytically that, individually, the Regge equations converge to zero 
as the second power of the lattice spacing, but that an average over local Regge equa¬ 
tions converges to zero as (at the very least) the third power of the lattice spacing. 
Numerical studies using analytic solutions to the Einstein equations show that these 
averages actually converge to zero as the fourth power of the lattice spacing. 


1 Introduction 

In 1961, Regge introduced a discrete form of General Relativity (which has since been named 
Regge calculus) 0] . Like General Relativity, Regge calculus is a geometrical theory of gravity. 
However, instead of using inhnitesimal distances (the metric) as the dynamical variables 
of the theory, Regge calculus employs finite distances as dynamical variables. Since its 
introduction, Regge calculus has been used for a variety of purposes (see [|]), but was 
initially envisioned as a tool for numerical relativity |l[]. Although work has been done in 
numerical relativity using Regge calculus, most applications share at least one of the two 
following undesirable (from the standpoint of practical numerical relativity) features: 

1. The method used is based on non-simplicial (non-rigid) blocks (i.e. prisms or hyper¬ 
cubes). 
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2. The method uses blocks with Euclidean signature. 

The hrst feature is undesirable for two reasons. First, not only must one specify extra 
“rigidifying” conditions on each block, but one must also make sure that the geometry on 
a face of a block matches up with the geometry of the face of an adjacent block (these 
“rigidifying” conditions have non-local constraints). Simplices, on the other hand, are rigid: 
hxing the values of the squared edge lengths completely determines the geometry inside 
the simplex, so there is no problem in trying to match up the geometry of two adjacent 
simplices. Second, it has recently been realized that Sorkin triangulations (triangulations in 
which the Regge equations decouple into single vertex evolutions) can be constructed from 
virtually any initial 3-dimensional simplicial manifold, thus forming a triangulation of the 
4-dimensional spacetime a ii- 

The second feature (Euclidean spacetime) is undesirable since the goal of numerical 
relativity is to model physically realistic processes. The fact is that once one has a Euclidean 
implementation of numerical Regge calculus, it is not an entirely trivial task to go over to 
Minkowskian signature.^ When taking derivatives of square roots of negative functions (as 
one does in Regge calculus with a —h ++ signature), one gets sign flips due to branch cuts 
in the square root function. One has to be very careful to avoid errors that are typified by 
the following “equation”: 

With a goal of using Regge calculus in numerical relativity, I have developed a code that 
calculates the Regge equations (and their derivatives with respect to the squared length of 
any link) on a simplicial manifold of arbitrary topology, with each simplex having a flat, 
Minkowskian metric. As a test for the code, I pose the following question: What is the 
order of the Regge calculus numerical method? That is, how fast do the Regge equations, 
when “evaluated” on a generic solution to the vacuum Einstein equations, approach zero as 
the lattice spacing approaches zero? (How one goes about “evaluating” the discrete Regge 
equations on a continuous solution of the Einstein equations is described in the next section.) 

In contrast to Regge calculus, the order of a given hnite differencing scheme is clear from 
the beginning: it is the order of the error introduced by truncating the Taylor expansion 
when approximating functions and their derivatives on a discrete lattice. For example, the 
derivative of a function /(x) defined on grid points Xj (evenly spaced; Xj+i — Xt = h) can be 
approximated at xq with a Taylor expansion about x = xq: 

fixi) = fixo + h) = /(xo) + fixo)h + 0{h^) 

^This is a technical point. Minkowskian signature introduces a bit of added complexity that must be 
dealt with, and it may not be as simple as adding a C-I--I- tag on timelike legs to an existing Euclidean code, 
as suggested in [|j. 
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namely, 


/-(,„). + 0(h) 

This is a first order approximation to the derivative of f{x) at xq. In a similar manner, any 
differential operator, £, acting on a set of functions, /(T), can be approximated by an 
order hnite difference operator, acting on functions / dehned on a lattice Xi characterized 
by uniform lattice spacing h. If /o(T) is a solution to the differential equation 

Cf = 0 

then, by construction, the hnite difference operator acting on fo{xi) (the continuum 
solution /o(x) evaluated at lattice sites Xi) will converge to zero as the power of h: 

4/o(x-;) = o + M" + o(h"+') 


where po is some constant. 

In the case of Regge calculus, it is not known, a priori, what order (in lattice spacing h) 
of error is introduced by discretizing the Einstein equations. This question is the subject of 
this paper. I will show analytically that, individually, the Regge equations, when evaluated 
on a continuum spacetime, converge to zero as the second power of the lattice spacing 
(there is a non-zero second order term).[| This result will be seen to be independent of 
whether or not the continuum metric satishes the Einstein equations. However, it will be 
shown analytically that averages (with weightings given by variations of each component 
of the continuum metric) over local Regge equations have a vanishing second order term 
only when the continuum metric satishes the Einstein equations. Studying the behavior 
of the third order term must be done numerically, using analytic solutions to the Einstein 
equations. These results show that, for the specihc solutions studied, the third order term 
of the average of the Regge equations also vanishes. 

2 Regge Calculus 

Regge calculus is based on the idea of approximating a curved continuum spacetime with a 
collection of hat simplices. For instance, a curved 2-dimensional surface can be approximated 
by a collection of 2-simplices (triangles). A general property of simplicial manifolds is that 
the curvature always resides on subspaces of codimension two. In the 2-dimensionaI case, 
this means that the curvature is concentrated on the vertices. In fact, the integral of the 
scalar curvature over a region of the simplicial manifold is simply a sum of the dehcit angles 
associated with each vertex. 

^This fact has been independently observed by Leo Brewin[^ 
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In four dimensions, the curvature is concentrated on 2-dimensional simplices (triangles). 
The Regge action is simply the Hilbert action evaluated on the simplicial manifold: 

S{1^) = J = ( 1 ) 

i 

where the sum on the right hand side is over all triangles (labeled by i), Ai is the area of 
triangle i, and is the dehcit angle associated with the triangle (boundary terms in 
the action are treated in |^). Here, we have set IQttG = 1. Demanding that the action 
be stationary with respect to small variations in the squared edge lengths yields the Regge 
equations. Each dynamical variable (squared edge length) has an associated Regge equation. 
The Regge equation associated with leg (i, j) is 

H cot = 0 (2) 

k 

where the sum is over all triangles {i,j, k) that have (q j) as one edge, Oijk is the angle Likj 
of triangle k), and eijk is the dehcit angle associated with triangle {i,j, k). 

One property of the Regge equations that will be important in the following sections is 
the fact that the equations are dimensionless. Only angles appear in the Regge equations. 
Specihcally, look at the expression for cotOijk in terms of the squared edge lengths 
and 1‘ji.: 

/2 /2 _ /2 

cot a ( 3 ) 

^/2i|4 + 211 , 1 % + 2111 % - /4 _ 4 _ 

We see explicitly that both the numerator and denominator have units of L^. In section ^ 
I will want to substitute a power series (in smallness parameter 6) for each R which has the 
form 

+ ■ ■ ■ ( 4 ) 

Observe that cotOijk (and, therefore, the Regge equation) would be indeterminate (0/0) 
when 5 = 0. In doing the actual calculations, equation (^) is scaled by a factor of 1/5^. The 
hrst non-zero term in a typical series expansion for P will now be of order unity, and the 

Regge equations will be perfectly well dehned at 5 = 0. Due to the complexity of the Regge 

equations, MATHEMATICA was used to perform the analytic calculations in section (^). 

3 A Formalism for Studying the Convergence Proper¬ 
ties of the Regge Equations 

In the past, various methods have been used to study the relationship between Einstein’s 
General Relativity in the continuum and the discrete Regge calculus. Both Sorkin and 
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Friedberg and Lee derive the Regge action from the Hilbert action by considering the 
limit of a sequence of continuum surfaces that approximate a piecewise flat surface. Barrett 
showed the equivalence between the Regge equations and the vanishing of the flow of 
energy-momentum across a certain hypersurface. Brewin |TI[ has shown the equivalence of 


the Regge action and the Hilbert action on almost flat simplicial spacetimes by treating the 
Riemann tensor as a distribution. 

The basic idea I use for comparing the continuum Einstein equations to the Regge equa¬ 
tions is the same as used in a paper by Cheeger et al [|^ (this idea can also be found in 
P]): assign the (signed) square of the geodesic length between two points in the continuum 
manifold whose metric satishes the vacuum Einstein equations to the squared edge length 

m 


of the corresponding link in the simplicial manifold. In 


it was shown (for Euclidean 


metrics) that the Regge action converges to the Hilbert action, in the sense of measures, 
as long as certain fatness conditions are satished by the simplices. Here, I ask a different 
question: how fast do the Regge equations converge to zero as the lattice spacing goes to 
zero? 


3.1 Evaluating the Regge Equations on a Continuum Spacetime 

For simplicity, this paper will use a hypercubic simplicial topology (hypercubes split into 
simplices). Consider a 4-manifold. At, equipped with a metric {g^i) whose signature is 
everywhere (—h ++). Furthermore, let be a solution to the vacuum Einstein equations 
Rfiu = 0. Let be coordinates on a region of At containing some point O with coordinates 
Xq = (0, 0, 0, 0). Now, introduce the hypercubic simplicial triangulation in the following way. 
First, construct a rectangular lattice in the coordinates {x^}. The lattice site (where 
n = (n°,n^,n^,n^)) is located at coordinates xh = (n° Ax°,n^ Ax^,n^ Ax^,n^ Ax^), where 
Ax^ are 4 arbitrary lattice spacing constants of order unity. To construct a triangulation 
from this rectangular lattice, we construct 15 links at each lattice site x^j: 


diagonal links rectangular links 




^n-\-hQ+h2) 




^n+ni+ri2) 

(^71: ^71+710 ) 

^n+ni+na) 

(^71: ^71+711 ) 

(Xfi^ 

(^717 ^71+712 ) 

^ri+no+iT 1+7^2) 

(^715 ^77+713 ) 

_^fii ^n+no+TLi+ris) 


^71+710+712+^3) 


_^fii ^71+711+712+713) 



(Xji, Xfi-\-fiQ+h\+h2+h^ 
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In this way, each hypercube is split into 24 4-siniphces. We now have a simplicial 
manifold with each vertex xa associated with a point on M. (the point with coordi¬ 
nates (n° Aa;°,n^ Ax‘^,n^ Ax^)). In order to study the limit of small lat¬ 

tice spacings, we multiply the coordinates corresponding to each lattice site with a 
“smallness parameter”, S. The coordinates associated with lattice site Xft are now 
(5 Aa;°, 5 Ax^, 5 Ax^, 5 Ax^). 

We now assign the (signed) square of the geodesic length of separation between points 
on the continuum manifold to the squared lengths of the links in the simplicial manifold 
(see Figure 1). For instance, to hnd the squared length of the link (a;(o,o,o,o), a;(i,o,o,o)) in the 
simplicial manifold. 


1. hnd the solution to the geodesic equation. 


(Px^^{X) 

dX^ 


+ Kpix{X)) 


dx°‘{X) dx^{X) 
dX dX 


0 


with boundary conditions a;^(A = 0) = (0,0, 0,0) and a;^(A = 1) = (5 Aa;°, 0,0,0), 
where A is an affine parameter ranging for 0 to 1. 


2. integrate the length of the geodesic : 


/ 



dx^^ dx^ 


' 9 


dX dX 


3. 


set the squared length of the link equal to (this quantity will be positive for spacelike 
geodesics, and negative for timelike geodesics). 


4 Analytic Results 

The geodesic equation cannot be solved analytically for an arbitrary metric. However, since 
all we want to do is hnd the hrst non-vanishing order of the Regge equations, we will be 
content with a power expansion (in S) for the squared geodesic distance of curves in the 
continuum manifold. 


4.1 Series Expansion for the Squared Geodesic Distances 

In general, we will need to hnd the power series expansion in 6 for the squared geodesic 
distances between two points Vi and Vf with respective coordinates {5 xf} and {5 x^} in 
terms of the metric at O, the Christohel symbols at O, the Riemann tensor at O, the hrst 
derivatives of the Riemann tensor at O, etc. Let x^{X) be a parametrization of the geodesic 
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with affine parameter A G [0, <5] (so that the affine parameter is of the same order as the 
length of the geodesic) that satishes the boundary conditions 


a;^(A = 0) = 6 Xi 
x^(\ = 5) = 5 X^j 


Now, expand x^{\) about the point A = 0: 


(5) 


= ,.(0) + _(0) A + ^^(0) A= + ^^(0) A3 +... (6) 

Taking this expansion, along with an expansion of the Christoffel symbols about O and 
substituting into the geodesic equation allows one to solve for the coefficients of eq. (§) in 
terms of the boundary values (the coordinates of the endpoints of the geodesic) and the 
Christoffel symbols and their derivatives at point O. One then expands the metric as a 
power series about O, and substitutes it into the hnal expression. 




dX 


dx^^ dx'" 


1 


dX dX 


/ dx^ dx^ A 

b-dA dA),./ 


(7) 


along with the power series of ^ (obtained by rearranging eq. (P)) to get the hnal power 
series expansion for the square of the geodesic distance between points Vi and Vf. 


4.2 Expanding the Regge Equations as a Power Series 

With a power series expansion for the squared geodesic distances in hand, we now evaluate 
the Regge equations on the simplicial manifold, where each squared length in the simplicial 
manifold is assigned the squared geodesic distance between the corresponding two points in 
the continuum manifold. A typical Regge equation can then be expressed as a power series 
in 6, 

'Rij{5) = 7^p(0) + '2\ ^ - 

where the coefficients of the series are functions of the metric at point (9, the Christoffel 
symbols at point O, the hrst derivatives of the Christoffel symbols at point O, etc. Recall 
that, due to the dimensionlessness of the Regge equations, in order to calculate the order 
coefficient in the expansion of a particular Regge equation, all of the expressions for the 
squared geodesic distances must be accurate up to and including the {n + 2)^^ order in 5. 
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4.2.1 Zeroth Order Results 


To calculate the zeroth order term in the expansion of a typical Regge equation, eq. (H), we 
require expressions for squared geodesic distances that are accurate up to and including the 


second order in 6. The analysis in section 4.1 gives us the answer: 


Ih = (9. 




- x’;)6^ + 0{6^) 


(9) 


Substituting these expressions into a typical Regge equation, and setting 5 = 0 is equivalent 
to evaluating the Regge equations on a continuum manifold with constant metric 
Therefore, all of the dehcit angles vanish, and the zeroth order term in the expansion of a 
typical Regge equation, eq. (||), is zero. 


4.2.2 First Order Results 


To calculate the hrst order term in the expansion of a typical Regge equation, eq. (||), we 
require expressions for squared geodesic distances that are accurate up to and including the 
third order in 5. Again, the analysis in section ^4.1| gives us the answer: 


+ ^/)(^/ - 4 )(^/ - ( 10 ) 


The hrst order term in the expansion of a typical Regge equation, has been com¬ 

puted on a number of different lattices (hypercubic lattice, randomly generated lattice, quan¬ 
tity production lattice [^), all with the same result: 


( 


dTZij 


) 


5=0 


0 


( 11 ) 


This result can be interpreted geometrically as a statement of the fact that there is no hrst 
order change in a vector that is parallel transported around a small closed loop. From the 
standpoint of numerical relativity, this result means that, at the very least, Regge calculus 
is a second order numerical method (Of course, the result is independent of whether or not 
the metric satishes the Einstein equations, since eq. (|I0|) depends only on the metric and its 
hrst derivatives (Christohels) at point O). 


4.2.3 Second Order Results 


Unfortunately, the problem of expressing the squared geodesic distances (in an arbitrary 
coordinate system), accurate up to and including the fourth order in 6 (needed to study the 
second order behavior of the Regge equations), in terms of the Riemann tensor (as well as 
the metric and Christohel symbols) at point O seems to be intractable. To remedy this, I 
will use Riemann normal coordinates |T^ about point O. One property of Riemann normal 









coordinates is that the Christoffel symbols vanish at point O. From eq. ([T0|), it can be seen 
that this causes the third order term in the expression for the squared geodesic distance to 
vanish. Using other properties of Riemann normal coordinates, namely, 



(12) 


(13) 

{dadpg^i/')Q T Rg^iav)Q 

(14) 


one can obtain the expression for the squared geodesic distance that is accurate up to and 
including the fourth order in 5\ 


l-f = 


2 {Rafifdu 






X- 


+ (15) 


One would hope that when substituting the above expression for squared geodesic dis¬ 
tances into the Regge equation expansion (^, one would observe the vanishing of the second 
order term when the Einstein equations are imposed = 0). Although there are cases 

(particularly when there is high lattice symmetry) where the second order term does vanish 
when the Einstein equations are imposed, it does not vanish in general. 


4.3 Averaging Over a Region of Spacetime 


From the above analysis, we see that only the zeroth and hrst order terms in the expansion of 
a typical Regge equation vanish. However, this result holds even if the continuum metric does 
not satisfy the Einstein equations!! Therefore, how conhdent can one be that Regge calculus 
is a good approximation to General Relativity? To motivate an answer to this, consider 
the analogous problem of approximating a continuous curve with a series of straight line 
segments. If we look closely at any one particular line segment, it does not look as if we are 
approximating the continuum curve very well. However, if we “step back” and look at many 
line segments, the approximation is better (see Figure 2). 

Therefore, we might expect a suitable average of the Regge equations contained in some 
region of spacetime to have a vanishing second order term when the Einstein equations are 
imposed on the continuum metric . How should one go about constructing this average? 


What should one use as a reasonable weighting? The continuum action principle suggests 
a partial answer [p!5|]. If we view the squared leg lengths in the Regge action as functions of 


the continuum metric. 




(16) 
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then the variation of the action with respect to variations in the metric {5gp.v) is simply 


ij \ 



(17) 


The hrst factor in the above sum, dS/dl^j, is the Regge equation that corresponds to leg 
(qj). The second factor can be used as the weighting of that particular Regge equation. 
There will be 10 independent weightings, corresponding to the 10 independent components 
of the metric. 

With the lattice constructed in section one could average over all Regge equations 
contained in a unit hypercube (it turns out that there are 65 of them in this particular case). 
I dehne TZavg as 

= (18) 

where the sum is over all legs (ij) contained in a unit hypercube, and rriij is the combinatorial 
factor 


rriij = 


number of 4-simplices inside unit hypercube that have (ij) as a leg 
total number of 4-simphces in the simplicial manifold that have (ij) as a leg 


(19) 


The 1/5^ in the weighting is to keep the weighting factor of order unity. There are ten 
independent weightings, each one corresponding to an independent component of the metric. 
Plugging in the hfth order accurate expressions for squared geodesic distances (^) into the 
Regge equations and expanding TZavgi^) as a power series about <5 = 0 reveals that the second 
order term of TZavgi^) vanishes when the Einstein equations are imposed at O. 


,(Pn. 


d6^ 


avg 1, _ g 

< 5=0 


( 20 ) 


This result holds for each of the ten independent weightings that correspond to variations 
of the ten independent metric components. 

One might wonder whether or not this result is an artifact of the usage of Riemann normal 
coordinates. To check, one can make an arbitrary coordinate transformation on equation 
(p!5|) before substituting the /^’s into the Regge equations. It turns out that the result, eq. 
(^), is independent of the coordinate transformation, thus, the result does not depend on 
the usage of Riemann normal coordinates. 


5 Numerical Results 

In this section, I study the convergence properties of the Regge equations numerically. In¬ 
stead of Ending a power expansion for squared geodesic distances (as in the previous section). 
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I numerically solve the geodesic equation for different analytic solutions of the vacuum Ein¬ 
stein equations (using the relaxation methods found in [0). I then evaluate the Regge 
equations on simplicial manifolds using these values. The first non-zero order term in a 
power expansion of or lZavg{5) can be found as the slope of the line of a ln(Jl{6)) 

vs. ln{l/6) plot. In all of the plots, the error bars are actually smaller than the data points 
appearing on the graph. This fact is a result of the accuracy of the relaxation methods used 
to solve the geodesic equations. I have verihed this by increasing the resolution of the grid 
points by more than one order of magnitude, and observing no change (to one part in 10®) 
in the values for the computed geodesic distances. 

Instead of picking on one Regge equation to evaluate and plot, I actually plot the sum 
of the absolute values of the Regge equations associated with the links in a unit hypercube 
(summing over the same links as in equation (|T8|)). I define TZabs to be 


TZabs = \TZi 


( 21 ) 




I also plot the weighted average of the Regge equations TZavg dehned by equation (|T^). 

The results from section ^ predict that TZabs should converge like 6^. Also, in section || 
it was shown that the zeroth. Erst, and second order terms of TZavg vanish. This leaves on 
open question: What is the first non-vanishing term of TZavg{S)Z 


5.1 Schwarzschild Metric 

The metric of the Schwarzschild solution has the form 

ds^ = - (^1 _ df + (^1 - + r^de^ + r^sin^e (22) 

I set M = 1, and put the origin of the lattice just outside the event horizon: Xq = 
(to, f'o, tto, 4’o) = (0, 3,7r/2, vr). The coordinate lattice spacing is /S.x^ = (At, Ar, A6*, A0) = 
(h, (5/2, ^/4, 5/4). The plot of TZahs is shown in Figure 3. The plot of TZavg, with a weighting 
induced by the variation of the time-time component of the continuum metric (the other 
weightings give the same results), is shown in Figure 4. 

The hgures show us a surprise. First, we would expect TZabs to converge as the second 
power in 5, but Figure 3 shows a fourth order convergence. Also, note the fourth order 
convergence of TZavg- To check these results, I did the same calculation, using Kruskal 
coordinates, and moving O inside the event horizon. In this case, the convergence properties 
are shown to be exactly the same as in Figures 3 and 4. 

Something to notice about both figures is the fact that, even though the lengths of 
the geodesics are of the order of the length scale defined by the curvature of the problem, 
the points corresponding to 5 = 1 are already on the converging line. This suggests that 
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Regge calculus may be superior to finite differencing in tracking the geometry of black hole 
spacetimes. 


5.2 Kerr Metric 


The Kerr metric, which is the vacuum solution corresponding to a rotating black hole, has 
the form 


ds'^ 


+ A‘^cos ‘^6 — 2 Mr\ 2 

r2 + A^cos^e J ^ 


AMArsin^O 


^2 _|_ \ 


-rj.v± j: XI oii± \ 7 7 / / ' I J- C/WO w \ 

'2 + A‘^cos‘^ 6 ) yr2 + ^2 — 2 Mr J 


\ / \ / \ 

/2 a 2 2n\jn2 / (r^ + A^f - (r^ + A^ - 2Mr)A^sm^d\ . 2 /, .,2 

f (r^ + A\os^e)de^ + ( ^- ,2 ^ 242cos 2^ -^-j 


(23) 


where M is the mass of the black hole and A is the angular momentum of the black hole 
divided by the mass: A = J/M. I set M = 1, A = 1/2 (so that the spacetime is asymptoti¬ 
cally predictable), and put the origin of the lattice at Xq = {toAo,(^o,<Po) = (0, 5,7r/2, tt). 
The coordinate lattice spacing is Ax^ = (At, Ar, Ad, Ac/) = (1.5(5, d, 0.2(5, 0.3(5). The plot of 
TZabs is shown in Figure 5, while the plot of TZavg, with a weighting induced by the variation 
of the time-time component of the continuum metric (the other weightings give the same 
results), is shown in Figure 6. 

As can be seen, this example shows TZabs converging as the second power of S, while TZavg 
still converges as the fourth power of S. 


5.3 Kasner Metric 

The Kasner metric, which corresponds to a spatially homogeneous anisotropic universe, has 
the form 

ds2 _ ^ ^2pi^^2 ^ ^ 2 p^^y 2 ^ ^2p3^^2 ^24) 

where the constants pi satisfy pi + P 2 + Ps = ^ and Pi + P 2 + P 3 = ^- I set pi = P 2 = 2/3, 
and ps = —1/3, and put the origin of the lattice at Xq = {to,xo,yo, ^o) = (2, 0,0,0). 
The coordinate lattice spacing is Ax^ = {At, Ax, Ay, Az) = {S,S,S,S). The plot of TZabs is 
shown in Figure 7, while the plot of TZavg , with a weighting induced by the variation of the 
time-time component of the continuum metric (the other weightings give similar results), is 
shown in Figure 8. 

We see that TZabs converges like the second power of 6 and TZavg converges like the fourth 
power of (5. These are the same convergence properties found with the Kerr metric. 
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6 Conclusions 


Both analytically and numerically, we see that the Regge equations individually converge as 
the second power of the lattice spacing, although this result is independent of whether or not 
the continuum metric satisfies the Einstein equations. We see that the second order term of 
an average over local Regge equations vanishes only when the continuum metric satishes the 
Einstein equations. Numerical studies (Figures 4, 6, and 8) show that, for certain analytic 
solutions of the Einstein equations, the third order term of IZavg also vanishes. Thus, Regge 
calculus can be used as a fourth order method in numerical relativity. 
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Figure Captions 


1. Method of “evaluating” a simplicial spacetime on a continuum space- 
time. 

2. Approximating a continuous curve with line segments. The approxima¬ 
tion is better if one looks at many line segments, rather than focusing 
in on any one particular segment. 

3. A plot of XniTZahs) vs. ln(l/h) for the Schwarzschild metric. 

4. A plot of XniTZavg) vs. ln(l/h) for the Schwarzschild metric. The 
particular weighting in TZavg (eq. 18) is the time-time variation of the 
metric. The graphs for the other 9 independent weightings are identical. 

5. A plot of XniTZahs) vs. ln(l/h) for the Kerr metric. 

6. A plot of XniTZavg) vs. ln(l/h) for the Kerr metric. The particular 
weighting in TZavg (eq. 18) is the time-time variation of the metric. 
The graphs for the other 9 independent weightings are identical. 

7. A plot of XniTZahs) vs. ln(l/h) for the Kasner metric. 

8. A plot of XniTZavg) vs. ln(l/h) for the Kasner metric. The particular 
weighting in TZavg (eq. 18) is the time-time variation of the metric. The 
graphs for the other 9 independent weightings have the same conver¬ 
gence rates. 
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